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Abstract 



Properties of spin polarized isospin asymmetric nuclear matter are stud- 
ied within the framework of the Brueckner-Hartree-Fock formalism. The 
single-particle potentials of neutrons and protons with spin up and down are 



o 

m 
o . 

, determined for several values of the neutron and proton spin polarizations 

o 

and the asymmetry parameter. It is found an almost linear and symmetric 
variation of the single-particle potentials as increasing these parameters. An 
pH _ analytic parametrization of the total energy per particle as a function of the 



asymmetry and spin polarizations is constructed. This parametrization is em- 



H ' ployed to compute the magnetic susceptibility of nuclear matter for several 

values of the asymmetry from neutron to symmetric matter. The results show 
no indication of a ferromagnetic transition at any density for any asymmetry 
of nuclear matter. 

PACS numbers: 26.60.-hc, 21.60.Jz, 26.50.-hx 

I. INTRODUCTION 

Pulsars are believed to be rapidly rotating neutron stars endowed with strong magnetic 
fields ||l|J^. Within the dipole magnetic model for pulsars the intensity of the surface (dipole) 
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magnetic field Bp can be inferred from the measured rotational period P of the pulsar about 
its rotation axis and from its time derivative P = dP/dt. For "canonical" pulsars [| this 
method gives values for Bp in the range 10^^-10^^ G. Despite the great theoretical effort, 
there is no general consensus regarding the mechanism to generate such a strong magnetic 
field in a neutron star. The field could be a fossil remnant from the one of the progenitor 
star. In fact, assuming magnetic flux conservation during the birth of the neutron star, a 
magnetic field of ~ 10^^ G could originate from the collapse of a main sequence star with a 
typical surface magnetic field of 10-10^ G. Alternatively, the field could be generated after 
the formation of the neutron star by some long-lived electric currents flowing in the highly 
conductive neutron star material. 

There are strong theoretical and observational arguments which indicate a decay of 
the magnetic field during the "life" of a neutron star P]. The physical processes which 
are responsible for the magnetic field evolution in neutron stars are not well understood, 
however, the observational data for pulsar distribution on the P-P plane can be reproduced 
by population synthesis codes (see e.g., ref. ^ and references therein quoted) in which the 
stellar magnetic field evolves according to the following exponential decay law 

B{t) = B^ + {Bo - 5oo)e-*/^^ , (1) 

where tb ~ 10^-10^ yr is the field decay time and -Boo ~ 10® G is a residue magnetic 
field. If the residue field is permanent and is not generated by some dynamo mechanism, 
it could be produced by a spontaneous ferromagnetic transition in the dense stellar core. 
Several authors have studied the possible existence of a ferromagnetic core in the liquid 
interior of neutron stars. First models, in which neutron star matter was approximated 
by pure neutron matter, were proposed just after the discovery of pulsars. Brownell and 



Here we use the term "canonical" to distinguish these pulsars from the group of "millisecond" 
pulsars, for which Bp ~ 10^-10^ G, and from a possible new family of pulsars (the so-called 
"magnetars") having Bp ~ 10^^ G. 
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Callaway and Rice |^ considered a hard sphere gas model and showed that neutron 
matter becomes ferromagnetic at kp ^ 2.3 fm~^. Silverstein ^ and 0stgaard found that 
the inclusion of long range attraction significantly increased the ferromagnetic transition 
density (e.g., 0stgaard predicted the transition to occur at kp ^ 4:.l fm~^ using a simple 
central potential with hard core only for the singlet spin states). Clark and Pearson 
and Saunier calculated the magnetic susceptibility for low densities {kp < 2 fm~^) using 



more realistic interactions. Pandharipande et al. ||T2|, using the Reid soft-core potential, 
performed a variational calculation arriving to the conclusion that such a transition was not 
to be expected for fci;' < 5 fm~^. Early calculations of the magnetic susceptibility within 



the Brueckner theory were performed by Backmann and Kallman employing the Reid 



soft-core potential, and results from a correlated basis function calculation were obtained by 
Jackson et al. with the Reid Vq interaction. A different point of view was followed by 



Vidaurre et al. , who employed neutron- neutron effective interactions of the Skyrme type, 
finding the ferromagnetic transition at /cj? ~ 1.73 — 1.97 fm~^. Marcos et al. have also 
studied the spin stability of dense neutron matter within the relativistic Dirac-Hartree-Fock 
approximation with an effective nucleon-meson Lagrangian, predicting the ferromagnetic 
transition at several times nuclear matter saturation density. On the contrary, no sign of such 



a transition has been found by Vidana et al. [|17l, who have recently studied the properties of 



spin polarized neutron matter within the Brueckner-Hartree-Fock approximation employing 
the realistic Nijmegen II and Reid93 interactions. In connection with the problem of neutrino 
diffusion in dense matter, Fantoni et al. |18| have recently employed the so-called auxiliary 
field diffusion Monte Carlo method (AFDMC) using realistic interactions (based upon the 
Argonne fis two-body potential [jl9[ plus Urbana IX three-body potential pO[), and have 
found that the magnetic susceptibility of neutron matter shows a strong reduction of about 
a factor of 3 with respect its Fermi gas value. They pointed out that such a reduction may 
have strong effects on the mean free path of a neutrino in dense matter and, therefore, it 
should be taken into account in the studies of supernovae and proto-neutron stars. 

In all these studies neutron star matter was approximated by a pure neutron gas. Never- 



theless, neutron star matter is also composed of protons, electrons, muons and other exotic 



constituents pl] , p2| . The importance of the presence of a small number of protons in neu- 
tron matter for the spin stability was pointed out by Kutschera and Wojcik [^]. These 
authors found that the ferromagnetic state, corresponding to completely polarized protons 
and weakly polarized neutrons, was energetically preferred over the nonpolarized one. Cal- 
culations performed by Bernardos et al. for strongly asymmetric nuclear matter within 
the relativistic Dirac-Hartree-Fock approximation confirmed also that the presence of an 
admixture of protons favors the ferromagnetic instability of dense matter. 

In this work we study the bulk and single-particle properties of spin polarized isospin 
asymmetric nuclear matter. Our calculations are based on the Brueckner-Hartree-Fock 
(BHF) approximation of the well known Brueckner-Bethe-Goldstone (BBG) theory of nu- 
clear matter. To describe the bare nucleon-nucleon interaction we make use of the nucleon- 
nucleon part of the recent realistic baryon-baryon interaction (model NSC97e) constructed 
by the Nijmegen group [^. We study the dependence of the single-particle potentials 
and the total energy per particle on the asymmetry parameter and neutron and proton 
spin polarizations. Further, we calculate the magnetic susceptibility for several values of 
the asymmetry parameter and in particular we explore the possibility of a ferromagnetic 
transition in the high density region relevant for neutron stars. 

The paper is organized in the following way. The formalism of the BHF approximation 
is briefly reviewed at the beginning of Sec. 0, whereas the calculation of the magnetic 
susceptibility is presented at the end of it. Section |T| is devoted to the presentation and 
discussion of the results obtained for the single-particle potentials, the total energy per 
particle and the magnetic susceptibility. Finally, a short summary and the main conclusions 
of this work are drawn in Sec. |V[ 
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II. FORMALISM 



Spin polarized isospin asymmetric nuclear matter is an infinite nuclear system composed 
of four different fermionic components: neutrons with spin up and spin down having densities 
PnT; and pni respectively, and protons with spin up and spin down with densities pp| and ppi- 
The total densities for neutrons protons (pp), and nucleons (p) are given by: 

Pn = PnT + Pni , Pp = PpT + Ppl ^ P ^ Pn + Pp- (2) 

The isospin asymmetry of the system can be expressed by the asymmetry parameter (3 — 
{Pn — Pp)/ Pi while the degree of spin polarization of the system is characterized by the 
neutron and proton spin polarizations Sn and S'p, defined as 

_ PnT - Pni _ PpT ~ Ppi (o\ 

i^n — ) i^p — ) K'^) 

Pn Pp 

Note that the value 5'„ = 5'p = corresponds to nonpolarized (i.e., = pn^ and pp^ = pp]) 
matter, whereas Sn — ±l('S'p = ±1) means that neutrons (protons) are totally polarized, 
i.e., all the neutron (proton) spins are aligned along the same direction. 

The single component densities are related to the total density and to the isospin and 
spin asymmetry parameters /3, and Sp via the equations: 

Pni = o Pn = 7;-P ' (4) 



Pni = ^ Pn = ^ ^P , (5) 

PpT = -^—Pp = ^P ' (6) 

1 — Sp 1 — Spl — (3 , , 

Ppi = -^—P^ = ^P ■ (^) 

Finally the Fermi momentum k'^ is related to the corresponding partial density by k'^ — 
(Gtt^Pi-ct)^^^, with T = n,p and a =T, i- 



Our calculations are based on the BHF approximation of the BBG theory extended to 
the case in which it is assumed that nuclear matter is arbitrarily asymmetric both in the 
isospin and spin degrees of freedom (i.e., pn-\ 7^ Pnj. 7^ PpT 7^ Ppi)- Therefore, our many- 
body scheme starts by constructing all the nucleon-nucleon G matrices which describe in an 
effective way the interaction between two nucleons (nn, np, pn and pp) for each one of the 
spin combinations (TT) Ti) it and [[). The G matrices can be obtained by solving the well 
known Bethe-Goldstone equation 

(fciTiai; A;2r2(T2|G'(u;)|fc3T3CT3; A;4T4CT4) = {kiTi(jr-,k2T2<J2\v\k3T3(j3; kiTiai) 
+ ^^(fciTicri; k2T2(T2\v\kiTiai; kjTjCTj) — -^{kmai; kjTjaj\G{uj)\k3T3a3; kinai) , 



where and am indicate respectively the isospin and spin projections (ritli),^?!!)) of 

— * 

the two nucleons in the initial, intermediate and final states, km are their respective linear 
momenta, v is the bare nucleon-nucleon interaction, Q^.^^^ .^.^^^ is the Pauli operator which 
allows only intermediate states compatible with the Pauli principle, and uj is the so-called 
starting energy. Note that the G matrices are obtained from a coupled channel equation. 
In practice a partial wave decomposition of the Bethe-Goldstone equation is performed, the 
Pauli operator and the energy denominator is replaced by an angle-averaged one, and the 
G matrices are solved using relative and center-of-mass coordinates. 

The single-particle energy of a nucleon with momentum k and spin projection a =T(i) 
is given by 

Era^^ + Uraik) , (9) 

where the single-particle potential Ura{k) represents the mean field "felt" by the nucleon 
due to its interaction with the other nucleons of the system. In the BHF approximation 

UTa{k) is calculated through the "on-energy shell" G matrix and is given by 

Uraik) = XI H Urar'a'ik) 

T'=n,po-'=T,i ^^Q^ 

= E II E Re{kTa;k'T'a'\G(u;^Era + Er'a')\kTa;k'T'a')A, 

r'=n,p o-'=T,i /c'<fej,'<^' 



where Urar'a'ik) is the contribution to Uraik) due to the Fermi sea of nucleons rV. Note 
that a sum over the Fermi seas of neutrons and protons with spin up and down is performed 
and that the matrix elements are properly antisymmetrized. We note also here that the 
continuous prescription has been adopted for the single particle potential when solving the 
Bethe-Goldstone equation. As shown by the authors of refs. P^ , |27[ ], the contribution to the 
energy per particle from three body clusters are diminished in this prescription. 

Once a self-consistent solution of Eqs. (||) and (0) is obtained, the total energy per 
particle is easily calculated: 

E 1 „ „ „ (n'k^ 1 ^ \ T V 

A = A^ ^ ^ [2^^2^^''^^n^A^A- ^^^^ 

This quantity is a function of p„f, p„|, and or, equivalently, of p, /5, and Sn and S^. 

The magnetic susceptibility x of system characterizes the response of this system to a 
magnetic field and gives a measure of the energy required to produce a net spin alignment 
in the direction of the field. In the case of nuclear matter it is defined by the 2x2 matrix 



1 

X 



np 



^ ^/Xpn '^/Xpp J 

where the matrix elements 1/Xij ^i^re given by 

1 dK 



(12) 



Here Aij is the magnetization of the species j per unit volume 

Mj = fijipj^ - pji) = fijPjSj , (14) 

with pj the magnetic moment of the species j, and Hi is the magnetic field induced by the 
magnetization of the species i, which can be obtained from 

= '^M- ■ ('^) 

Using Eqs. ([l^ and (|l^) the matrix elements l/xij can be written as 
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1 

Xij 



P 



'd\E/A)\ 



(16) 



fJ'iPifJ'jPj \ dSidSj J 

where the second derivatives can be taken at Si = Sj = if the magnetic field is assumed 
to be small. 

It is convenient to study the magnetic susceptibility of the system in terms of the ratio 
det(l/x)/det(l/xF)5 where xf is the magnetic susceptibility of the two component free 
Fermi gas. Writing the total energy per particle as a sum of the kinetic, T/A, and the 
potential, V/A, energy contributions, the ratio between both determinants can be written 

as 



det(l/x) 



1 



d'^jV/A) 
d'^{T/A) 



/ 



1 + 



V 



82 (T/A) 



f d^V/A) 



d'^{T/A) 



(17) 



det(l/xF) 

The stability of matter against spin fiuctuations is guaranteed if det(l/x)/det(l/xF)> 0, 



indicating a change of sign of the ratio the onset of a ferromagnetic phase in the system. 



III. RESULTS 

A. Single- particle potentials 

In Fig. |l| we show the neutron Uni (upper panels), and proton f/p|, f/pj (lower 

panels) single-particle potentials evaluated at p = 0.17 fm~'^ for several values of Sn, Sp, and 
f3. Results for symmetric matter {(3 = 0) are plotted on the left and middle panels, whereas 
results for asymmetric matter with f3 = 0.5 are reported on the right ones. Solid lines show as 
a reference the results for nonpolarized matter {Sn = Sp = 0), while dotted lines refer to spin 
polarized matter. On the left panels we show the effects on Uraik) of a partial polarization 
{Sn = 0.75) of the neutron component in proton- unpolarized {Sp = 0) symmetric matter. 
As it can be seen, the single-particle potentials of both species split off when a partial 
polarization of the neutron spin is assumed. The single-particle potential for neutrons with 
spin up (down) becomes less (more) attractive with respect to Un{k) for nonpolarized nuclear 
matter. For the proton single-particle potential, we find an opposite qualitative trend being 

8 



f/pl (Upi) more (less) attractive than the proton single-particle potential for nonpolarized 
matter. The additional partial proton polarization (middle panels) produces on Uraik) 
a reverse global effect with respect to the one produced by polarizing neutrons. This is 
particularly evident in the case of Upa{k) (compare the results in the lower left and lower 
central panels). Finally, to obtain the results presented in the right panels of Fig. we 
further introduce an isospin asymmetry {f3 = 0.5) in the system. Notice that the isospin 
asymmetry generates a splitting of the neutron and proton single particle potentials even in 



the case of spin-unpolarized matter p8 . 

The splitting, Ur^-Uri, in the neutron and proton single-particle potentials can be mainly 
ascribed to two reasons: (i) the change in the number of pairs which the nucleon under 
consideration \k, r, a) can form with the remaining nucleons \k' < kp'^' , r', a') of the system 
as nuclear matter is spin polarized, and (ii) the spin dependence (and isospin dependence for 
asymmetric spin polarized matter) of the nucleon-nucleon G matrices in the spin polarized 
nuclear medium (see Eq. (||)). 

In the general case of spin polarized asymmetric matter, the effect of the neutron and 
proton spin polarizations and the isospin asymmetry on the single particle potentials can be 
clarified by considering their partial contributions UraT'a'ik) in such a way as to single out 
explicitly the dependence on their respective phase space: 

Un-] = Un^n^ + Un}n[ + t^nTpT + ^"Tpi = Qn^MPn'i + Qn^nlPnl + 5'nTpTPpT + dn^plPpl 5 (1^) 

Unl = Unln-\ + f^jni + f^nipj + ^"ipi = 9nln^pn'\ + 5'ninJ.Pni + 5'riMPpT + dnlplPpl , (19) 

Up\ = Up^n] + Upi^ni + ^pTpT + ^pTpi = 9p\n\Pn\ + dp^niPni + fi'pTptPpT + 9p\piPpi , (20) 

Upi = UpM + Uplnl + Uplp-^ + Uplpl = gpinlPnl + QplnlPnl + fl'pipTPpT + 9plplPpl , (21) 

where Qrar'a' is the average value of the matrix element {kra; k'r' a'\G\kTa] fc'r'o"')^ in the 
Fermi sphere with radius k' < kp"'' and the density factor Pt-Zq-/ in each term arises from the 
integral over the corresponding Fermi sea. 
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For small values of the asymmetry parameter << 1) and the spin polarizations 
{\Sn\, \Sp\ « 1) one can neglect the dependence on /5, Sn and Sp of the average G matrices 

{k, p) and 

9nW ~ 9nlnl ~ fi'pTpT ~ 9plpi = dl , 

9n'(nl ~ gnln-] ~ Qp^pl ~ Qplp^ = 92 , ^^^^ 
9n]p-^ ~ 9nlpl ~ 9p-\nt ~ 5'pini = 93 , 
9n]pl ~ 5'nipT ~ 9pUl ~ 5'pinT = 94 ■ 

Clearly, even under these assumptions, in the most general case 9i 7^ 92 93 9i because 
they receive contributions from different spin and isospin channels. Whereas gi receives 
contributions only from the spin and isospin triplet {S = 1,T = 1) channels, g2 has in 
addition contributions from the spin singlet ones, from channels with 5 = 1 and T = 0, 1, 
and g4, from channels with S* = 0, 1 and T = 0, 1. 

Using Eq. (|22|) the single-particle potentials can then be rewritten in terms of the average 
G matrices grar'a' as: 



Un^ ~ ^ [{gi + g2 + g3 + g^) + {gi + 5-2 - 5-3 - gA)l3 

+ {g, - g2)il + (3)Sn + (^?3 - gA){l - (3)Sp] , 



Uni ~ ^ [(5-1 + g2 + g3 + gi) + {gi + 5-2 - 5-3 - 6-4)^ 
-((71 - g2){l + /3)5„, - ((73 - g4){l - P)Sp] , 



^pt ~ ^ [(^1 + ^2 + ^3 + ^4) - (^1 + g2-g3- gA)P 
+ig3 - ^74)(l + P)Sn + (gi - g2)il - P)Sp] , 



Upi ~ ^ [(^1 + ^2 + ^3 + ^4) - {gi + g2-g3- gi)l3 



(23) 



(24) 



(25) 



(26) 



-{g3 - ^?4)(i + P)Sn - (gi - g2)ii - l3)Sp] , 

where Eqs. (^), (|^) (^ and (|^ have been used to write pra in terms of p, /?, S'„ and Sp. 
These equations show explicitly the dependence of the single-particle potentials on the spin 
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polarizations and isospin asymmetry. This dependence is tested in Fig. where the value at 
= of the single-particle potentials Un-\, Uni, Up^ and Upi at p = 0.17 fm~'^ is plotted as a 
function of neutron (left panels) and proton (right panels) spin polarizations for two values 
of the asymmetry parameter: j3 = (upper panels) and (3 = 0.5 (lower panels). In order to 
make more clear the discussion Sn{Sp) is taken equal to on the right (left) panels. The 
above equations predict a linear and symmetric variation of the single-particle potentials on 
/3, Sn and Sp. This prediction is well confirmed from the microscopic results reported on 
Fig. ^, although deviations from this behavior are found at higher values of the asymmetry 
and spin polarizations. These deviations have to be associated to the dependence on j3, Sn 
and Sp of the average G matrices grar'a' which has been neglected in the present analysis 
(see Eq. (0)). From the above expressions it can be seen that under charge exchange (i.e., 
doing the changes P — /?, Sn — ^ Sp, Sp Sn) the role of neutrons and protons with the 
same spin projection is exactly interchanged. This can be clearly seen by comparing the 
left and right top panels of the figure, which correspond to the particular case of symmetric 
matter. It is also clear that neutrons (protons) with spin up and down interchange their 
roles when a global flip of the spins is performed (i.e., changing Sn by —Sn and Sp by —Sp). 



B. Energy per particle 

The total energy per particle of neutron, asymmetric, and symmetric matter is shown 
on the left, middle and right panels of Fig. ^, respectively. In the three panels solid lines 
show results for nonpolarized matter, whereas those for totally polarized matter are re- 
ported by dotted lines. Note that in totally polarized asymmetric and symmetric matter, 
we have distinguished two possible orientations of the neutron and proton spins: that in 
which neutron and proton spins are aligned along the same direction (i.e., Sn = Sp = ±1), 
and that in which neutron and proton spins are orientated along opposite directions (i.e., 
Sn = ±1, Sp = =f1)- As can be seen from the figure, totally polarized matter is always more 
repulsive than nonpolarized matter in all the density range explored for any value of the 
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asymmetry parameter. Note also that the case in which neutron and proton spins are paral- 
lelly orientated is less repulsive than the case in which they have an antiparallel orientation. 
To highlight the effects of the nuclear interaction on the variation of E/A as nuclear matter 
is polarized, we plot in Fig. ^ the kinetic (upper panels) and potential (lower panels) energy 
per particle for nonpolarized and totally polarized matter. The kinetic energy contribution 
of totally polarized neutron matter, asymmetric or symmetric matter is always larger than 
the corresponding one of nonpolarized matter, simply due to the fermionic nature of nucle- 
ons. Also, the potential energy contribution is always more repulsive in the totally polarized 
case. This can be understood by considering separately the contribution to the potential 
energy per particle of the spin singlet and spin triplet channels. In Tables | and |I| we show 
these contributions for nonpolarized and totally polarized neutron and symmetric matter 
at densities p = 0.17 fm~'^ and p = 0.4 fm~'^, respectively. Note, firstly, that an impor- 
tant amount of binding is lost in totally polarized neutron matter and in totally polarized 
symmetric matter with neutron and proton spins parallelly orientated due to the absence in 
these cases of the contribution of the spin singlet channels. Note also, in the totally polarized 
symmetric case, that when the orientation of neutron and proton spins is antiparallel this 
contribution, although it is present, is much less attractive or even repulsive (see Table ^ 
than the corresponding one of nonpolarized symmetric matter. In addition, in all cases the 
contribution from spin triplet channels in totally polarized matter is always less attractive 
(or even repulsive at high density) than the corresponding one to nonpolarized matter. In 
particular, it is mainly this contribution the one which explains the difference between the 
energies of totally polarized symmetric matter with a parallel and an antiparallel orientation 
of neutron and proton spins, being in the antiparallel case less attractive or more repulsive. 
An interesting conclusion which can be inferred from our microscopic calculations is that a 
spontaneous transition to a spin polarized state, i.e., to a so-called ferromagnetic state, of 
nuclear matter is not to be expected for all the possible isospin asymmetries ranging from 
symmetric to pure neutron matter. If such a transition would exist, a crossing of the energies 
of the totally polarized and nonpolarized cases would be observed in neutron, asymmetric 
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or symmetric matter at some density, indicating that the ground state of the system would 
be ferromagnetic from that density on. As can be seen in Fig. there is no sign of such a 
crossing and, on the contrary, a spin polarized state becomes less favorable as the density 
increases. 

It would be interesting to characterize in a simple analytic form the dependence of the 
total energy per particle on the asymmetry and spin polarizations. The kinetic energy 
contribution is already analytic and well known. It is given by 



jiP, P, S^, S,) = \^\ [(1 + Pf'i^ + Snf" + (1 + Pf'il - S^fl' ^^^^ 

+ (1 - (3fl\l + Spfl^ + (1 - f3f/\l - Spf/^] , 

being kp = (37r^p/2)^/^ the Fermi momentum of nonpolarized isospin symmetric matter. An 
idea of the possible terms appearing in the potential energy contribution can be extracted 
from our previous phase space analysis of the single-particle potentials. We start from the 
BHF approximation of the potential energy per particle V/A defined through Eq. (|TT|) and 
making use of Eqs. ([18|)-(^TD for the single particle potential Uraik), we can write 

T-yrEE E Wa'(^,P)PrV' • (28) 
^ T.o- r',(T' k<ky^ 

Next we use Eq. (|2^) and perform the average of the quantities gi{k, p) over the correspond- 
ing Fermi sphere with radius ky . Finally, we arrive, after some little algebra, to 

^[^1 - -g,]{l + PfSl + ^[^1 - ^2](1 - PfSl + ^[^3 - ^4](1 - P')SnS, , 

where cji are just the averages values of the quantities gi. Following this simple analysis we 
can finally infer the form of the total energy per particle 

f (p, /3, 5„, S,) = ^(p, /3, 5„, S,) + K)(p) + ^i(p)/3' + ^2(p)(l + PfSl ^ ^ 

^ ^ (30) 

+V2{p)il - fifSl + V,{p){l - P^)SnS, . 

This parametrization is consistent with the spin and isospin structure of the nucleon- 
nucleon interaction, in the sense that if we consider a particular configuration, (3, Sn, Sp, of 
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the system, a global flip of the spins does not change the energy, whereas a flip only of all 
neutron or protons spins does. It is also clear that due to charge symmetry, a system with 
a conflguration, (3' = —(3, S'^ = Sp, Sp = Sn, will have the same energy as the original one 
(note that the coefficients of the terms (1 + P^S^ and (1 — P)^Sp are the same). 

The coefficients Vo{p), Vi(p), V2(p) and V^i^p), whose density dependence is shown in Fig. 
^, have been determined in the following way 

Vo{p) = j{p,P = 0,Sr, = 0, Sp = 0) , (31) 
Viip) = j{p,P = l,Sn = 0, Sp = 0)- Voip) , (32) 
V^ip) = ^{p,(3 = 0,Sn = 1, Sp = 0)- Voip) , (33) 

Vsip) = j{p,P = 0,5„ = l,Sp = l)- Voip) - 21^2(p) . (34) 

It is clear, however, that their determination is not unique. Choosing them in this way, we get 
a parametrization which reproduces with a good quality (see Figs. |] and |^ and the discussion 
below) the results of the BHF calculation of the total energy per particle for values of /?, Sn 
and Sp around their values for nonpolarized symmetric matter (i.e., (3 = 0, Sn = 0, Sp = 0). 
Note that for nonpolarized matter this parametrization reduces to 

jip, P, 0, 0) = jip, P, 0, 0) + Voip) + Viip)P' , (35) 

which can be identified with the usual parabolic approach of the nuclear matter energy per 
particle if the kinetic energy contribution is expanded up to order P^. 

In order to test the quality of this parametrization, we show in Fig. ^ the total energy per 
particle at p = 0.8 fm~^ as a function of the proton spin polarization for different values of the 
neutron spin polarization and two values of the asymmetry parameter: P = (right panel) 
and P = 0.25 (left panel). Circles, squares, diamonds and triangles show the results obtained 
from the BHF calculation, whereas those obtained from the parametrization are reported by 
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solid lines. As can be seen from the figure, the dependence on the spin polarizations and the 
asymmetry parameter predicted by Eq. ( pO| ) is well confirmed from the microscopic results. 
It is interesting to note that for a fixed value of [3 and Sn (with Sn 0), the minimum of the 
energy happens for a value of Sp ^ 0. However, this is not an indication of a phase transition 
to a ferromagnetic state, because the real ground state of the system, as we have seen, is that 
of nonpolarized matter in all the range of densities and isospin asymmetries considered. For 
completeness, we compare in Fig. ^ the results for the total energy per particle as a function 
of density for three arbitrarily asymmetric and spin polarized situations: (3 = 0.25, Sn = 
0.3, Sp = 0.4, p = 0.5, Sn = 0.5, Sp = 0.25, and p = 0.8, Sn = 0.6, Sp = 0.2, obtained from 
the BHF calculation and from the parametrization. As in the previous figure, symbols show 
results obtained from the BHF calculation, while sold lines refer to those obtained with the 
parametrization. The quality of the parametrization is quite good, as can be seen in both 
figures, with deviations from the microscopic calculation of at most 7 — 8% only for the 
combinations of /5, Sn and 5*^ with the largest values. Higher deviations are found, however, 
at higher values of these parameters, being this an indication that the coefficients need to 
be refitted when considering matter with |/?| ~ 1, l^nj ~ 1 and l^pl ~ 1. 



C. Magnetic susceptibility 

The ratio det(l/x)/det(l/xi?) can be evaluated in a very simple analytic way from Eq. 
(p!7|) if the parametrization of Eq. (|30|) is assumed, giving 
det(l/x) , 12m1/2(p)r,^ , mi/3 , /i mi/3i , ^em^ 

(36) 



with m the average mass of the nucleons. This ratio is shown in Fig. ^ as a function of the 
density for several values of the asymmetry parameter /3 from neutron to symmetric matter. 
Note that an important increase of the ratio happens as soon as a small fraction of protons 
is introduced in the system. This can be understood by looking at Eq. (|36D and Fig. ^ from 
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which it is clear that the ratio det(l/x)/det(l/xF) increases as the asymmetry parameter 
13 decreases from 1 to 0. According to the criteria for the appearance of a ferromagnetic 
instabihty, such an instabihty should appear when det(l/x)/det(l/xF)= 0. Nevertheless, it 
can be seen from the figure that det(l/x)/det(l/xF) increases monotonously with density, 
and a decrease is not to be expected even at higher densities for any value of the asymmetry 
parameter. Therefore, it can be inferred from these results that there is no sign at any 
density of a possible ferromagnetic phase transition for any asymmetry of nuclear matter. 

As we said previously, the coefficients of the parametrization should be refitted when 
large values of the asymmetry parameter and the spin polarizations are considered. In 
particular, for the case of pure neutron matter it is more convenient to redefine coefficient 
V2(p) in the following way 

V2{P) = \ QU/^ = = = 0) - ^(p,/5 = l,^n = 0,5p = 0)) , (37) 

in order to get a better value of the magnetic susceptibility. We show in the figure by 
open circles the resulting neutron magnetic susceptibility obtained by refitting coefficient 
V2(p) according to Eq. (|37|) . We find a good agreement with the results of ref. jl^ where 
this coefficient is fitted in this way. For comparison, we show also the results for neutron 
matter obtained recently by Fantoni et al. (filled circles). We note here, that our results 
confirm the reduction of about a factor 3 of the magnetic susceptibility of neutron matter 
with respect to its Fermi gas value found by these authors, being the differences ascribed 
to the method (AFDMC) and the potential (AU6' two-body + UIX three-body) employed 
by them. It seems, therefore, that this reduction is largely independent of the two-body 
modern potential used, and, in addition, that the effects of the three-nucleon interaction on 
the magnetic susceptibility are not large. 



IV. SUMMARY AND CONCLUSIONS 



Employing a realistic modern nucleon-nucleon interaction (NSC97e) we have studied 
properties of spin polarized isospin asymmetric nuclear matter within the Brueckner- 
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Hartree-Fock approximation. We have determined the single particle potentials of neutrons 
and protons with spin up and down for several values of the neutron and proton spin po- 
larizations and the asymmetry parameter. We have found that the potentials exhibit an 
almost linear and symmetric variation as a function of these parameters. Deviations from 
this behaviour occur at higher values of the asymmetry parameter and spin polarizations. 
These deviations have to be attributed to the dependence of the nucleon-nucleon G matrices 
on /3, Sn and Sp. 

We have calculated the total energy per particle as a function of the density for totally 
polarized and nonpolarized neutron, asymmetric and symmetric matter. We have found 
that in the range of densities explored (up to 7po) totally polarized matter is always more 
repulsive than nonpolarized matter for any asymmetry. This is due to a combination effect 
of the kinetic and potential energy contributions. We have also found in the totally polarized 
case that asymmetric and symmetric matter is more repulsive when neutron and proton spins 
are antiparallelly orientated than when all the spins are aligned along the same direction. 

We have constructed an analytic parametrization of the total energy per particle as 
a function of the asymmetry parameter and spin polarizations. The quality of this 
parametrization has been tested, finding deviations from the microscopic calculation of at 
most 7 — 8%. Nevertheless, deviations are higher for higher values of the asymmetry param- 
eter and spin polarization, and the coefficients of the parametrization need to be refitted 
when large values of the spin and isospin asymmetries are considered. 

Employing this parametrization we have determined the magnetic susceptibility of nu- 
clear matter for several values of the asymmetry from neutron to symmetric matter in terms 
of the ratio det(l/x)/det(l/xi;'). We have found that this quantity increases monotonously 
with density, from which it can be inferred that a phase transition to a ferromagnetic state 
is not to be expected in nuclear matter at any density for any asymmetry. 

Finally, our results confirm the reduction of about a factor 3 of the magnetic susceptibility 
of neutron matter with respect to its Fermi gas value found recently by Fantoni et al. |T^]. 
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TABLES 

TABLE L Contribution of the spin singlet and triplet channels to the potential energy per 
particle of neutron and symmetric matter at p = 0.17 fm~^. Second and third columns show results 
for nonpolarized and totally polarized neutron matter, respectively, whereas those for nonpolarized 
and totally polarized symmetric matter are reported on the fourth, fifth and sixth ones. Partial 
waves have been included up to total angular momentum J = 4. 





Neutron matter 




Symmetric matter 






5„ = 


Sn = ±1 


Sn — Sp — 


Sn — Sp = zbl Sn 


= ±l,5p = Tl 


5 = 


-21.970 




-14.414 




-0.616 


5 = 1 


-0.247 


9.960 


-24.650 


-23.173 


-7.647 


Total 


-22.217 


9.960 


-39.064 


-23.173 


-8.263 



TABLE IL As in Table | but for p = 0.4 fm^^. 





Neutron matter 




Symmetric matter 






5„ = 


Sn = ±1 


Sn — Sp — 


Sn — Sp = zbl Sn 


= ±l,5p = Tl 


5 = 


-37.222 




-22.619 




8.291 


5= 1 


10.090 


48.173 


-33.102 


-20.014 


3.413 


Total 


-27.132 


48.173 


-55.721 


-20.014 


11.704 
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FIGURES 

FIG. 1. Neutron (n|,n|) and proton (p], pi) single-particle potentials as a function of the 
linear momentum A; at p = 0.17 fm^"^ for several values of /?, Sn and Sp. Left and middle panels 
show results for symmetric matter, whereas those for asymmetric matter are reported on the right 
one. Solid lines show as a refence the results for the nonpolarized case, while dotted ones refer to 
the spin polarized cases. 

FIG. 2. Neutron (n| and nj.) and proton (p| and pi) single-particle potentials at A; = and 
p = 0.17 fm^'^ as a function of neutron (left panels) and proton (right panels) spin polarizations. 
Upper (lower) panels show results for symmetric (asymmetric) mattter. Sn (Sp) is taken equal to 
on the right (left) panels. 

FIG. 3. Total energy per particle as a function of the density for nonpolarized (solid lines) and 
totally polarized (dotted lines) neutron (left panel), asymmetric (middle panel) and symmetric 
(right panel) matter. 

FIG. 4. Kinetic (upper panels) and potential (lower panels) energy contributions to the total 
energy per particle as a function of the density for nonpolarized polarized (solid lines) and totally 
polarized (dotted lines) neutron (left panel), asymmetric (middle panel) and symmetric (right 
panel) matter. 

FIG. 5. Density dependence of the coefficients of the parametrization defined in Eq. (pO|). 

FIG. 6. Total energy per particle at p = 0.8 fm~'^ as a function of the proton spin polarization 
for different values of the neutron spin polarization and two values of the asymmetry parameter: 
P = (right panel) and (3 = 0.25 (right panel). Circles, squares, diamonds and triangles show the 
results of the BHF calculation, whereas solid lines refer to the parametrization defined in Eq. (|30|). 
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FIG. 7. Total energy per particle as a function of the density for three arbitrarily spin po- 
larized asymmetric situations: /? = 0.25, Sn = 0-3, Sp = 0.4, (3 = 0.5, Sn = 0.5, Sp = 0.25, and 
(3 = 0.8, Sn = 0.6, Sp = 0.2. As in Fig. ^, symbols show the BHF results, whereas solid lines 
correspond the parametrization defined in Eq. (|30|). 

FIG. 8. Ratio between the determinats of the matrices 1/x and 1/xf as a function of the 
density for several values of the asymmetry parameter f3 from neutron to symmetric matter. 
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